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ABSTRACT 

Combined  Computational  Fluid  Dynamics 
Reynolds  Averaged  Navier-Stokes  and  LES 
model  calculations  were  completed  for 
comparisons  with  wind  tunnel  measurements  of 
flow  along  the  afterbody  and  in  the  base  region  of 
a  circular  cylinder  aligned  to  a  Mach  2.5 
freestream.  This  cylinder  is  representative  of  a 
missile  at  zero  degrees  angle-of-attack. 
Comparison  data  includes  afterbody  and  base 
pressure,  high  quality  non-intrusive  Laser 
Doppler  Velocimetry  and  Pressure  Sensitive  Paint 
measurements,  and  turbulent  kinetic  energy. 

1.  INTRODUCTION 

The  flow  in  the  base/separated  region  of  a 
tactical  missile  is  very  complex  even  though  the 
geometry  is  quite  simple.  There  are  several 
reasons  this  is  the  case.  First,  the  external  flow 
over  the  missile  body  is  supersonic  while  the  flow 
in  the  base  region  is  mixed  supersonic/subsonic 
flow.  Second,  this  flow  is  characterized  by  shock 
wave/boundary  layer  interactions  at  the  lip  of  the 
base  region.  Third,  the  flow  is  turbulent  in  the 
base/separated  region  and  highly  transient  in 
nature. 

There  has  been  a  long  history  of  attempts  to 
accurately  describe  this  flow  going  back  to  the 
1950’s  where  integral  models  were  used  and 
continuing  to  the  2000’s  where  sophisticated 
computation  fluid  dynamic  models  were  used. 
However,  none  of  these  models  have  produced 
results  that  provide  accurate  descriptions  of  the 
flowfield. 

Due  to  the  technological  advances  that  have 
been  made  in  recent  years  in  the  field  of  Large 
Eddy  Simulation  (LES),  parallel  processing 
methods,  and  high  speed  computers,  it  was 
determined  that  it  was  time  to  make  another 
attempt  to  solve  these  elusive  problems. 


Quality  data  sets  extending  from  the  early 
measurements  of  Reid  and  Hastings1  through  the 
JANNAF  series2,3  and  the  more  recent 
measurements  by  Dutton  and  students4"11  have 
consistently  challenged  the  computational 
capability  of  the  fluid  dynamics  community  and 
even  the  recent  calculations  by  Papp  and  Ghia12 
show  disappointing  comparisons  with  the 
cylindrical  blunt  base  data  of  Herrin  and  Dutton6. 
Additional  work  using  more  advanced  turbulence 
models  (EASM,  LES)  by  Ayyalasomayajula  13, 14 
also  failed  to  show  improvement  in  the 
comparisons. 

The  objective  of  this  work,  then,  is  to 
determine  if  recent  improvements  to  the 
Hybrid/RANS/LES  AMRDEC  CFD  flow  field 
model  calculations  show  any  better  agreement 
with  the  UIUC  LDV  measurements.  Of  interest  is 
the  visualization  of  the  flow  field  about  a 
cylindrical  body  at  0  degrees  angle-of-attack, 
surface  pressure  along  the  after  body  and  base 
surfaces,  and  the  flow  field  properties  in  the  base 
region.  A  significant  addition  to  these 
comparisons  is  the  inclusion  of  the  entire  wind 
tunnel  facility  in  the  calculational  domain  since 
small-scale  facilities  can  affect  the  results  due  to 
flow  interactions  with  the  facility  geometry. 

2.  EXPERIMENTAL  FACILITY 

The  experimental  data  described  herein  were 
collected  in  the  Gas  Dynamics  Laboratory  at  the 
University  of  Illinois  at  Urbana-Champaign.  This 
axi- symmetric  wind  tunnel  facility  was 
specifically  constructed  for  the  investigation  of 
base  flow  phenomena.  A  photograph  of  the  wind 
tunnel  is  shown  in  Figure  1.  High  pressure  dry  air, 
as  illustrated  in  Figure  2,  enters  the  stagnation 
chamber  from  a  regulated  supply,  passes  through 
a  honeycomb  flow- straightening  module,  and 
expands  through  a  converging-diverging  nozzle 
with  central  sting  to  an  inviscid  Mach  number  of 
2.5  (actual  measurements  indicate  a  freestream 
Mach  number  of  2.46).  This  uniform  supersonic 
flow  (to  within  1%)  then  passes  through  a  test 
section  and  supersonic  diffuser  to  exit  the  facility. 
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The  objective  of  these  tests  was  to  simulate 
base  flow  about  a  missile  at  zero  degree  angle-of- 
attack.  A  cylindrical  base  configuration  was 
fabricated  from  stainless  steel  and  was  mounted 
centrally  through  the  wind  tunnel  nozzle.  A 
schematic  of  the  afterbody  and  flowfield  is  shown 
in  Figure  2. 

The  test  section,  essentially  a  box  to  mount 
four  flat  optical  windows  for  optical  access,  forms 
very  complex  recirculating  cavity  flows,  which 
are  quite  difficult  to  model.  However,  according 
to  Dutton,  it  was  found  that  adjusting  the 
stagnation  pressure  such  that  the  nozzle  exit  lip 
pressure  matched  the  pressure  in  the  test  cell 
could  reduce  waves  that  emanate  from  the  nozzle 
exit  lip  which  disrupt  the  near-wake  flowfield. 
Under  these  conditions,  interference  waves,  as 
viewed  in  Schlieren  imagery  (fig.  3),  were  seen  to 
vanish  at  the  matched  pressure  condition.  Full 
details  of  the  wind  tunnel  and  model  construction, 
as  well  as  the  LDV  measurements,  are  given  in 
References  15  and  16. 

3.  COMPUTATION  METHODOLOGY 

The  AMRDEC  CFD  models  the  full  Navier- 
Stokes  (FNS)  equation  set  providing  an  aero- 
thermo-chemical  plume  /  airframe  predictions  for 
unsteady-flows17.  The  CFD  code  numerics 
include  1D/2D/Axi/3D  finite  volume 
discretization  with  an  implicit,  higher-order 
upwind  (Roe/TVD)  formulation  (2nd  order  for 
RANS  and  5th  order  for  LES).  The  turbulence 
model  includes  a  k-e  formulation  with 
compressibility/vortical  upgrades  over  standard 
models  and  includes  several  low  Re  near-wall 
formulations.  A  variable,  one/two  equation 
Prandtl  number/Schmidt  model  has  been  recently 
added.  All  equations  are  solved  either  fully 
implicitly  or  loosely  coupled. 

This  Reynolds  averaged  code  has  been 
extended  to  handle  three  types  of  LES  modeling. 
The  first  is  Miles/LES  that  has  no  subgrid  or  wall 
modeling  capability,  the  second  is  a  Wall 
Blended/LES  model  that  uses  a  hyperbolic 
blending  function  between  RANS  and  LES  near 
walls  and  the  third  model  is  a  point-by-point 
Hybrid  RANS/LES  model  (AJ).  This  AJ  model 
alters  from  pure  energy  transferring  (LES)  to 
energy  extracting  (RANS)  through  a  point-by¬ 
point  blending  function.  The  RANS  solution 
provides  a  coarse  limit  of  the  LES  solution.  This 
hybrid  model  estimates  the  local  range  of  scales 
before  applying  the  blending  via  two  additional 
equations.  The  first  equation  estimates  the  extent 


of  unresolved  flow  and  the  second  equation 
estimates  the  largest  length  scale  of  the  flow. 
These  two  equations  are  similar  in  form  to  the 
RANS  k-8  formulation  but  have  different 
meanings.  At  each  point  in  the  flow,  the  RANS 
and  LES  viscosity  is  computed.  A  blending 
function  that  is  based  on  the  cell’s  location  (i.e. 
wall  boundary  layer,  shear  layer,  free  stream,  or 
mixed  condition)  is  than  used  to  compute  the 
local  viscosity. 

Boundary  conditions  for  this  after  body/base 
flow  problem  were  specified  as  follows.  Uniform 
subsonic  conditions  were  fixed  at  the  stagnation 
chamber  inflow  boundary  for  total  conditions  of 
517  kPa  pressure  and  293  K.  Gradient 
extrapolation  procedures  were  used  at  the  outflow 
boundary.  Surface  wall  conditions  were  viscous 
with  imposed  adiabatic  no  slip  flow.  Boundary 
layers  were  grid  resolved  to  <  1  Y+  using  a  So- 
Zhang-Speziale  wall  model.  For  comparison, 
RANS  only  calculations  reset  both  the  turbulent 
kinetic  energy  and  production  of  turbulence 
values  to  chamber  background  values. 

Several  numerical  grids  were  developed  to 
take  advantage  of  the  symmetry  of  the  wind 
tunnel  configuration.  This  included  the  use  of  a 
pie  shaped  grid  defining  the  stagnation  chamber, 
throat,  nozzle,  and  center  sting.  This  pie  wedge  is 
then  imposed  onto  a  full  360  degree  multi-block 
(interface)  grid  just  upstream  of  the  test  section  to 
allow  for  any  possible  downstream  flow  feed¬ 
back.  The  remaining  grid  is  used  to  model  the 
test  section.  The  test  section  is  modeled  without 
walls  -  just  outflow  boundaries. 


4.  RESULTS  AND  DISCUSSION 

Three  different  grid  types  were  used  to 
compute  the  solutions  for  comparison  with  data. 
The  first  grid  included  a  full  360  degree  region 
from  chamber  to  test  section  and  had  2.6  million 
points  with  one  point  every  3  degree 
circumferential  angle.  The  second  grid  also 
included  a  full  360  degree  region  from  chamber 
to  test  section  and  had  9.2  million  points  with  one 
point  every  1.5  degree  circumferential  angle.  The 
third  grid  made  use  to  the  “pie  wedge”  grid 
mentioned  above  to  model  the  chamber  through 
nozzle  region.  This  grid  had  a  total  5.4  million 
points  but  due  to  the  large  grid  point  savings 
around  the  center  sting,  this  grid  had  the  most 
points  in  the  base  and  downstream  shear  layer. 
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Using  the  three  grids  with  the  three  different 
LES  models  produced  the  following  results: 

1)  For  base  pressure,  only  the  AJ  model  gave 
a  near  constant  base  pressure  profile  as  shown  in 
the  comparison  with  measurements  (fig.  4).  Grid 
size  was  more  important  than  LES  model 
selection  in  terms  of  base  pressure  correlation. 
The  best  base  pressure  comparison  results  were 
produced  using  the  AJ  model  run  with  5.4  million 
grid  points.  The  error  was  +6%.  All  other 
models  produced  results  where  base  pressures 
were  below  the  measured  pressures  and  ranged  in 
error  by  as  much  as  -25%.  The  Miles  case  was  the 
most  numerically  stable  run  (highest  stable  time 
step)  while  the  wall-blended  model  was  the  least 
stable. 

2)  Centerline  Turbulent  Kinetic  Energy  (TKE) 
predictions  varied  greatly  from  the  measured 
values  (fig.  5).  The  9.2  million  grid  case  gave  the 
best  results.  This  figure  shows  that  the  measured 
TKE  values  on  the  centerline  do  not  vary  greatly 
over  all  axial  locations  (5000  to  7500  m2/s2). 
However,  measured  off-axis  values  are  much 
higher. 

3)  For  the  centerline  axial  velocity,  all  models 
predicted  the  same  axial  profile  shape  except  for 
the  Miles  model  with  2.6  million  grid  points  case 
(fig.  6).  Most  models  predicted  the  zero  velocity 
point  downstream  of  measured  point  and  most 
models  predicted  a  larger  negative  axial  velocity 
component  in  the  recirculation  zone. 

4)  No  radial  fluctuations  were  evident  along 
the  centerline,  however  the  wall  blended  model 
predicted  some  noise  near  the  base  face.  The 
measured  data  recorded  a  variation  in  radial 
velocity  on  the  centerline  of  -10  to  +20  m/s  (fig. 
7). 

Axial/radial  velocity  and  TKE  measurements 
off  the  centerline  were  compared  in  detail  with 
the  AJ  model  using  5.4  million  grid  points  (figs. 
8-11).  The  axial  velocity  comparison  clearly 
shows  a  much  longer  base  flow  region  than 
measured.  The  radial  velocity  prediction 
compares  well  with  the  measurements.  The 
velocity  vector  plot  clearly  shows  the  predicted 
base  region  (as  defined  by  the  zero  velocity  line) 
as  larger  and  longer  than  measured.  The  TKE 
measurements  show  a  much  lower  value  over  the 
region  of  the  base  flow.  The  predicted  mixing 
layer  profile  clearly  does  not  match  the  measured 
shear  layer. 

Looking  to  improve  the  modeled  solution,  it 
was  thought  that  the  tunnel  walls  might  impact 
the  base  flow.  An  axisymmetric  RANS  case  was 


run  with  and  without  the  presence  of  the  tunnel 
walls  (for  this  case  the  windows  were  not 
modeled).  The  measured  base  pressure  ratio  was 
was  found  to  average  0.56  P/Pinf.  For  the  case 
where  the  walls  were  not  modeled,  the  RANS 
normalized  predicted  base  pressure  was  found  to 
be  0.41  P/Pjnf  (-27%  error).  For  the  case  where  the 
walls  were  modeled,  the  RANS  predicted  base 
pressure  was  found  to  be  0.54  P/Pinf  (-4%  error)  - 
a  significant  improvement.  Clearly,  the  presence 
of  the  tunnel  walls  affects  the  axisymmetric 
RANS  solution.  The  effects  of  the  tunnel  walls 
on  a  LES  prediction  were  not  addressed  in  this 
study. 

5.  SUMMARY  AND  CONCLUSIONS 

Inclusion  of  the  tunnel  geometry  improves  the 
agreement  with  data. 

1 .  The  AJ  RANS/LES  point-by-point 
modeling  gave  the  best  base  pressure 
results.  The  model,  however,  did  an 
inadequate  job  of  predicting  the  size  and 
extent  of  the  downstream  base  flow.  This 
model  needs  more  work. 

2.  The  gross  features  of  base  flow 
phenomena  are  correct  but  local  features 
differ  significantly  from  the 
experimental  measurements.  LES 
predictions  for  base  pressure  are 
sensitive  to  both  the  model  and  grid 
points. 

3.  RANS/LES  wall  blended  modeling  is 
presently  inadequate.  The  current 
hyperbolic  wall  blending  function 
between  full  RANS  and  full  LES  is  not 
physical. 

4.  Miles/LES  modeling  of  the  cylindrical 
base  is  inadequate.  It  is  not  meant  for 
and  should  not  be  used  for  base  flows 
where  boundary  layers  are  present. 

Overall: 

It  appears  that  the  model  does  not 
contain  sufficient  physics  to  handle  the 
shock  wave/boundary  layer  interactions  that 
are  predominant  in  the  base  flow  cases 
studied. 
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Figure  8.  TKE  Comparison  to  Measurements 
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Figure  9.  Axial  Velocity  Comparison  to  Measurements 
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Figure  10.  Radial  Velocity  Comparison  to  Measurements 
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Figure  11.  Vector  Velocity  Comparison  to  Measurements 
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